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Abstract. This paper considers the limiting distributions of 7r CT , the sta- 
tionary distribution of infinitely-many-alleles diffusion with symmetric selec- 
tion. In Tj the limiting distribution of was obtained except critical case, 
<t = A log 8, A = k(k + l),fc > 1. This article gives a positive answer to the 
critical problem. Therefore, all limiting distributions of n a with ovcrdominant 
selection and small mutation are obtained. 



1. Introduction 

The most extensively studied model in population genetics is neutral infinite- 
allele model, where mutations always generate completely new allele types. If we 
use x = (xi, X2, ■ ■ ■ ) to represent the allele frequency in decreasing order, then 
the infmitely-many-neutral-alleles diffusion [3] is the associated diffusion process 
characterized by generator 

i A , o 2 ^ a 



where = {(xi, x%, ■ ■ ■ )\x\ > X2 > • • ■ > 0, x i — !}• Poisson-Dirichlet 
distribution, hereafter denoted as PD(0), is its stationary distribution. If symmetric 
selection is considered, we will end up with the inhnitcly-many-alleles diffusion with 
symmetric selection introduced by Ethier and Kurtz in [4], characterized by 



J") Xi(Xi - ip2(x)) ■ 



G a = G + (7^1,(1, - <^ 2 (x)) — ,x e v c 

i=l 

where V?2(x) = Yl^-i x 1- F° r a given allele frequency x, the homozygosity, denoted 
by H2, is defined to be ^(x). Since the allele frequency is random, then H2 is a 
random variable as well. 

The infinitely-many- alleles diffusion with symmetric selection has stationary dis- 
tribution 7T CT defined as 

ir a (dx) = C a exp{aip2{x)}PD(0)(dx),x g V w , 

where C a is a normalized constant. Choosing symmetric selection is not only bio- 
logically important [12J . but also mathematically convenient for the corresponding 
stochastic process is Markovian. If a > 0, then the selection is of underdominant 
type; if a < 0, then the selection is overdominant. 

Mutations and selection are two important evolution forces in population genet- 
ics. It is commonly accepted that mutations contribute substantially to the genetic 
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variability. Overdominant selection, however, is also suggested to be the crucial fac- 
tor for maintaining biological diversity, please refer to and references therein. 
When mutations and overdominant selection are both large, the effect of overdom- 
inant selection is hardly pronounced; while when mutation is small and overdomi- 
nant selection is large, the effect of overdominant selection can be observable. In [9] 
Gillespie conjectured that overdominant selective model is indistinguishable from 
neutral model based his numerical results. This conjecture was certified by Joyce, 
Krone and Kurtz |10j , Dawson and Feng [T] . In [7J , large deviation principle (here- 
after LDP) of PD(9) with small mutation is obtained, the speed is log ^ and the 
LDP rate function is 

{0, x € Li 

n—1, x G L n ,x n > 0, n > 2 
+oo, x £ L 

where, L n = {(xi,--- ,x n , (),■■■) G V : Yn=i x i = 1 }> L = U^Li L n- Here J(x) 
exhibits some properties similar to energy ladder structure. As an application, the 
LDP of ir a was also considered. When a = Xlog8 the LDP rate function is 

S\(x) — J(x) + \ip2(x) — inf { h n — 1 : n > 1}, 

n>i n 

which is essentially different from that of PD(6>), thus, the effects of overdominant 
selection are pronounced. It was observed in [7) that when A € (k(k — 1), k(k + 
^ 1) the limiting distribution of 7r CT is S(X r .. i,o,— )! f° r the critical case 
A = fc(fc + 1), k > I, the LDP rate function S(x) has two zero points, but the law 
of large numbers remains open, for more details please refer to [5] . 

The main result of this paper provides a positive answer to the law of large 
numbers at critical value. When a — A log 9, A = fc(fc + l),fc > 1, the limiting 
distribution of ir a is obtained. It is 

oo 

5 m-i)Mk+i)] ( x ) s (h- .bo,-)- 

k=i 

Therefore, the limiting distribution of homozygosity under Tr a is also obtained and 
it is 

oo 

(!) <5(fc(fc-i),fc(fc+i)] (A)(5 1 . 

fe=i 

As a matter of fact, ((TJ actually rigorously confirms the simulation results in 
where the distribution of heterozygosity, defined to be H2 = 1 — H2, concentrates 
at 5 or according to the simulation results when mutation is small and ratio of 
selection to mutations is large; thus the distribution of homozygosity concentrates 
at \ or 1. But based on results in [7], if we consider the ratio of selection to the 
logarithm of mutations, then the simulation results may be better. 

Furthermore, the probability density function of homozygosity under PD(#) ex- 
hibits spikes at 4, k > 1, as well, though its limiting distribution concentrates only 
at 1 as 9 — > 0. This statistical property of homozygosity was observed by Stewart 
in [6] . In [5] Higgs notices that there are some similarities of frequency distributions 
between population genetics and statistical physics, and some simulations in [5] also 
show that homozygosity has spikes at > 1- Interestingly homozygosity turns 
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out to be an important quantity in those statistical physics models as well. In [2] 
those spikes of homozygosity distributions were shown to be the common features 
of three typical physics models: random map model, mean field spin glass model 
and random broken object problem. As Higgs argued in [5], however, homozygosity 
can hardly have physical realizations in those models; on the contrary, homozygos- 
ity is a directly measurable quantity in population genetics. Due to results in this 
article, the singularities of homozygosity can be even balanced out under mutations 
and ovcrdominant selection. Even though there are similarities between models in 
population genetics and physics, I am still not sure if we can regard mutations as 
functions of absolute temperature and treat overdominant selection as some exter- 
nal field, then A becomes the ratio of this field to the logarithm of temperature, 
then it will exhibit countable many phase transitions when temperature is near 
and the field is strong. 

As we know, critical cases in physics are notoriously difficult, therefore, it is 
not an easy job to figure out the limiting distributions of homozygosity when A = 
k(k + l),fc > 1. In this article, by calculating the asymptotic moment generating 
function of homozygosity under Ti a , the limiting distribution of homozygosity is 
obtained for all A > 0; in addition, by making use of LDP of ir a , the asymptotic 
distribution of n a is obtained as well. 

The whole paper is organized as follows. In section 2, we will present the main 
theorem on limiting distribution of homozygosity and its proof. In section 3, proofs 
of all lemmas will be shown in details. 

2. Main results 

Theorem 2.1. For a = X\og6, A > 0, the limiting distribution of homozygosity, 
H2, under the distribution ■n a is 



Remark 1. As you can see, the limiting distribution of homozygosity is pretty 
much a Dirac comb. At each critical point, the distribution stay the same as its 
previous one, after the critical point, the distribution jump to another one, such a 
phase transition phenomenon resembles the metastability. The idea of proof is quite 
simple, that is, we are going to figure out the limiting moment generating function of 
homozygosity under -k g . Thus, the limiting distribution of homozygosity is obtained. 

Before we present the proof, we need the following lemma, the proof of which 
are presented in the next section. 

Lemma 1. The moments of heterozygosity = EH\ = E(l — H2) k has the 
following form 




k 



m k = 
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where 

2k + 9 2 k klT(k + 9) 



Ak ^- 2k T(2k + l + 6) 
. ... ^ 2k + 6 2 k k\ T{k + l + 6) A 

A ^») = E —Yn r(2k + i + e) Al "-^ 



/=p-i 

Let's define Ak, p = Ak. p (0), then 

_ 2 fc fc!(fc- 1)! 
M ~ (2fc)l 

and 

^ 2 fc /c!r(fc + /) , 

Ak -= £ i f(2fc+TT ' p ~ 1 ' p " 2 - 

Then Afc iP does not depend on anymore, and we have the following lemma 
their upper bounds and relations. 

Lemma 2. J/ we p, for large k, V# G [0, 1], we have 

^A ktP <A ktP {e)<A k ,p<^, Vfc>p>l; 
|A fe)J ,(0) - A fciP | < 6»pA fcjP , I <p < k. 
Lemma 3. For A > 0, we have 



f e l ± i -^^A k , l (e) = o. 

e^o z — ' ^ k\ 

J=[A]+1 fe=i 



Lemma 4. Suppose that a n ,b n are two positive sequences, and lim„_ ) . 00 
then lim^+co ^= ob "^k = c. 
Lemma 5. For p > 1 fixed, we have 

a r 1 ( p \ k 

where C\ — y/w, and C p+ \ = Cp^/nC^-)^ . 

Lemma 6. Define C k j = ^Zo («) {^) l A k - s , h then 



XI . 1 I 



{X-l)(l + l) ; ki a i + r 
Lemma 7. For A > 2, define 

WA] nl^oo (Alogl)" . ... 

El^E^^l^W 

and 



Af/^ _ 2^=1 " 2^fc=Z fc! ^fc+n,i 

EI^E^^I^m : 

then when u{u — 1) < A < u{u +1), u > 2 7 we have 

lim = lim K*(0) = C^-) n . 
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[PROOF OF THEOREM El]: 

Proof. Let's denote (f>H to be the moment generating function of homozygosity H 2 
under Tt a . Thus, 

, > _ £exp{-tff 2 } ■ exp{A log 6H 2 } 
0H[t) - Eexp{X\og9H 2 } 

For some technical reason, we need to multiply the common term A log ^ to the 
numerator and denumerator, then 



<j>n(t) = e 



_ t Ee^{t{l - H 2 )} • exp{Alog i(l - H 2 )} 



-e 



£exp{Alog±(l-ff 2 )} 
, , , ^ f" E(l - H 2 y exp{A log |(1 - H 2 )} 



n=l 



111 



oo 



£exp{Alog|(l-tf 2 )} 

00 (A log i)'" _ „ \n+m 



t 2 = rm=0 m! ^H 1 ~ ^ 2 

n! v^°° (a log iy» jj v 



If we denote the limit of 4>ii{t) by ipnit), then we have 

00 n (A log |) fc /i fflW 

= e"* 1 + hm £ V Efc=1 ~grS ,=1 ; * + " ,,( ) ) 

i + EL^Eti^W 

By Lesbegue's dominant convergence theorem, we can switch the order of summa- 
tion and limit. Thus, 



1> H (t) = e-*(l + T ^(lim S 2^1=1 A ^,iW0 } 

^ w V ^ n! v e-*o x j_ y-°° (*i°g?) fc ^k 



' ' 1 — y — ( iijii — 

OO (Alog^) fc ^fc + n . (Q\Ql, 



Now we claim that 

„, 1 v^oo (Alog ^) fc v^ fc +™ a ta\ai 
Um m " + 2^ fc= i gj 2^=i+fc A k+n ,i {0)0 = 

1 + 2^fc=l k\ ZW=1 A k,l{0W 

Indeed, we have 



1 1 V°° (Alog|) fc y^fc A 
00 /at 1 \ k+n 

fe=l " ' l = l+k 



6 



YOUZHOU ZHOU 



By lemma [21 we have 

fc=i z=i+fc fc=i z=fc+i 

00 (A§log^ fe 



fe! v v 2 7 '2- 
fc=i 

' 3 \n % if'„. > '# lo, 



(1 _ (-)™)(e A 2 1o « t - 1) -> 0, as 6> -> 0. 



2 6> v v 2 

We have used the fact that m n — > 0, as 9 — > 0, which is due to PD{9)(dx) — > 
5(1^,...) (dec), as — > 0. By Lesbegue's dominant convergence theorem, we can 
switch the order of summation and limit, thus we have 



=e- t (l + V^lim^ 

n=l J- 



i=l# Efc=i fel ^fe+n,z(fj 



Now what we need to show is, for u(u — 1) < A < u(u + l),u> 1, we have 
(2) hm ^t^^-^A^M = u-1 

^i+e^e^t^mw ' " 

Once we have got the above equation, then 

^( t ) = e -'(l + f^(^)«) 
z — : n\ u 

n=l 



U 



i n! 

Thus, 4>H{t) — E^Li ${u(u-i),u(u+\)] " ■ Therefore, the limiting distribution 
of H 2 under 7r CT is Y^Li *(u(«-i),u(u+i)](A)5i. 

Now we are going to certify the claim ©. Firstly, when < A < 2, we have 

< hm ^M3C^i ^-W^) 
1 + 1^1=1 V l^k=l fel A k,l{°) 1=1 k=l 

fe=l ' fc=2 

1=3 k=l 

We can actually show that the above three terms approach as 9 — > 0. Indeed, by 
lemma [31 we have 

1=3 k=l ' 1=3 k=l 

< 4 y(-)'e A1 °^ = -i(-) 3 = -9 3 - x -> 0, as 0^ 0. 
- 6> A 2 2 

2=3 
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and by lemma [21 we have 



0S9 g(AMr Aw(9)s9 g(AMf Ai+j>i 

fc=l ' fc=l 

Y^oo (log 
Z^fc=o fc! 

The above limit is due to lemma [5] and lemma [4] Similarly, 

fc=2 ' fc=l 



< ^ M 2^ ^o, as 8^ 0. 

- y^oo (21og^) fc 

2^fc=0 fc! 



Thus, we have for < A < 2, 

lim S^'E^^^Wfl) = Q = i-i 
^i + E^E^t^mW 1 

Secondly, for u(u — 1) < A < u(it + 1), tt > 2, then A > 2. we can show that 
(3) lim ^ k = l f* +n '' W = lim AT* (0) , 

»-° i + E,=! & J2Zi ^V^Aw (*) e "° 

then by lemma [3 we have proved claim ([2]) . Now we only need to verify ([3]) . To 
this end, 

W-M fli V°° ( A1 °S?) fc 4 /m , V™ nl v^°° (A log ^) fc . , fl , 
_ l^i=i V 2-rfc=j fc! ■^fc+n,;(,fe | j + Z^=i+[a] » Z^fc=; g A fc+rt j(frj 

i+E~ifl'Er=i ^p^w ' 

then by lemma [31 we know 

Eoo nl v^oo (A log |) fc a fn\ 

,=i+ ^ gSfc=, ~r:f fc+mi,(fl) -.o, as^o. 

1+E^Efcl* ^t^MW 
And by lemma [2] and lemma El we have 

1=1 k=l ' fe=l 



Y^oo (Alogl) fc . 
l^k =l fc! ^fc+ n.l 

Y^oo (log 
Z^fc=0 fc! 



> — — i , fc > +oo, as -> 0. 



Similarly, E|=i EfcH; ^1 9 " ' (g) -> +oo, as -> 0. Therefore, claim 1 is 
proved. □ 
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Theorem 2.2. For a given a = A log 6, A > 0, as 9 — »■ 0,71V converges weakly to 
the following 

oo 

XI <W-l),fc(fc+l)] ( A ) 5 (£,- ,£,0,- )■ 
fe=l 

Proof: We only need to consider the critical case A = (k + l)k, k > 1. For a fixed 
fc > 1, and A = (k + l)k, we need to prove that 7r CT converges weakly to S ( i ... i 

For a given / e C(Voo), thcnVe > 0,36 < fc(fc+ \ )+1 , such that Vx e • • • , ±,0, 

we have 

|/(^)-/(p---,p0,---)l<e. 

Thus, 

| ^ /(^(^-/(I,... ,p0,---)| 
-I I (.f(x)-.f(i 0,--.)K(dx)| 

< / v l/W-/(J,-,pO,-)kW 

= / i/w-/(t,--- ,io,---)i^w 

Js x <s k k 

<2\\f\\ 0O 7r a (S x >5)+ f |/ (a; ) _/(!,. .. ,^,0,...)\^(dx) 

J(s x <s)n(\<p 2 -i\>6) k k 

+ / l/w-/^-,^,-)!^) 

J(SA<5)n(^2-i|<5) fe fe 

< 2||/|| 00 ( 7 r M (S A > 5) + n xfi (\V2 - \\ > 6)) 

+ [ \m-f(l, ■■■ ,p o,---)K(dx) 

J(S x <6)n(\ V2 -±\<6) fc fc 

By LDP of 7IV and the weak convergence of H2 under ir a , we have 

lim ir a (Sx >S) = lim n a (\<p 2 - T \ > 6) = 0. 
Moreover, we claim that 

(4) (S x <6)n(\cp 2 -±\<6)CB 5 (1,--- ,p0,-")- 

Then we have, 

/" 11 
limsupl / f(x)n a (dx) -/( 7 ,--- , -,0,---)| < 

Letting e — > 0, we have 

lim / f(x)Tr a {dx) =/(-,••• , -,0, •••)• 

Therefore, 71V converges weakly to 



c. 
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Now we need to show the claim (j3}, for A = (k + l)k, since 

S x (x) = J(x) + {k + l)kLp 2 {x) - inf{ ^ + 1 ^ + n - 1}, 

n>l n 

then 

Sx\l„{x) =n-l + k(k + l)ip 2 \ Ln - 2k. 
Because <P2\l„(x) has a unique minimum point (K ■ ■ ■ , i, 0, • • • ), then 

Sx\ Ln (x) >n-l + fc(fc + 1) -2fc. 

n 

By the monotonicity of the righthand function in n, we know it attains its minimum 
at k and k + 1, and for n ^ fc, fc + 1, 

Sa|l„(*) > min{fc - 1 - 1 + *£±i> - 2fc, fc + 2 - 1 + ^til _ 2 fc} = JL_ > 5. 
Then (5 A < 6) = (S\ < 6) n (L& U L fe+ i), thus, 



(5 A < 6) n - i| < 5) 
l 



fc 1 



(5 A < 8) n (|¥> 2 - r | < 5) n i fe u (5 A < <5) n (|<^ 2 - -| < 5) n L fc+ i 



i 



A- 



But 5 A |L fc+1 = *(* + l)(^|L fc+1 - ^i), then (5 A < 5) n L k+1 = (<p 2 < + 
^L) n L fe+ i. Since 5 < k{k + l)+l , then + ^ < ± - <J, thus we have 

Va; € Lfc+i, 

^ 2(x)< fc(fcTT) + fcTT < ^" (5< ^ (x) - 

Therefore, (5 A < 5) n (|<p 2 — || < 6) C\ L k+1 = 0. We have 

(5a < S) n (|^ 2 - i| < 5) = (5 A < s) n (|^ 2 - i| < 5) n L k 

Since Vx e L k , \ipi{x) — jM < 8, we have | — 6 < Yli=i x 1 < 7 + <^> an d 

_ < ! < 1- _. 

A/ A/ A^ Aj^ A; 



Therefore, y pr - f < mini<,-< A a^ < maxi<j< fc ^ < v/jr + f- Then Vx G 

(5 A <«J)n (|^2-^| <^)ni t , 

^,(p-,io,-)) = E J ^ 



fe' 

^ i i i i 

< > — max I a:< — — I = (1 r) max |a?i — — , 



i=l 



.. , 1 <5 1. . / 1 S l n 



= max{ 
< 6. 



8 5 
k k 



I ill /i i S i 1 

k 2 k ^ k V k 2 k k 
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Therefore (5a < 5)n(|y>2-|| < S) C B s (j,--- >i' >"'): claim ® is thus proved. 



3. Proof of Lemmas 

Now we will embark on a long journey to prove the previous lemmas, but we will 
omit the proof of lemma |4] for it can usually be found in a standard textbook. 
[PROOF OF LEMMA Q]: 

Proof. Dchnc Vi = U 1: Vi = (1 - U x ) ■ ■ ■ (1 - Ui-i)U u i > 2, where {U u i > 
1} are i.i.d. Beta(l,0). Then (Fi,V 2 , •■•) follows GEM distribution. Denote 
(V(i) , V(2) , • ' ' ) the order statistics of ( Vi , V2 , • • • ) by decreasing order, then ( V(i) , V( 2 ) , 
follows Poisson-Dirichlet distribution. Since H 2 = J2iLi^(i) = Eti^ 2 > then 
1 - H 2 = (1 - f/f) - (1 - Uxf + (1 - f7i) 2 (l - £T 2 ) = 2U X (1 - U x ) + (1 - (7i) 2 (1 - H 2 ), 
where H 2 — YnLi ^> an d Vi = ^2, Vi = (1 — ?7 2 ) •■•(! — E/i)[/i + i. We can see that 
(Vi, V2, • • • ) follows GEM as well and is independent of Ui. Thus, £7(1 - H 2 ) k = 
E(l - H 2 ) k , and 



m fe = £7(1 - H 2 ) k = {2Ux{\ - U x ) + (1 - t/i) 2 (l - H 2 )) k 
= Ej2( k l )(2U 1 (l-U 1 )) k - l (l-H 2 ) 1 



1=0 

( k l )E(2U 1 (l^U 1 )) k - l E(l-H 2 ) 1 

1=0 



E 



1 

' 'fc\ 2 fc -'r(fc - / + i)r(fc + 1 



1 ) r(2k + 1 



-mi. 



1=0 

If we isolate m^, we have 

_ Q l ^ 2k + 6 2 k k\ r(k + l + 8) 2k + 9 2 k k\T(k + 9) 

(5) mk ~ jri 2fc 2'/! r(2fc + l + 0) m ' + 2fc f(2jb+T+d) ' 

where fc > 2, and mi = j^g- We claim that has the following expansion 

k 

1=1 

where A k ,i = Z)«=i-i fsjg rfffc+T+e) Indeed, for fc = 1, this is obvious. 

Assume that mk-i has the above expression, then for rrik, by ([5]). we have 

^ 2fc + fl2 fc fc! r(fc + z + g) ' 2fc + e 2 fc fc!r(fc + g) 

^ 2fc 2'/! r(2fc+l + 0)^ iM ^ 2fc r(2fc + l + 0) 

1 ' 1 2fc + 9 2 fe fc! r(fc + I + 6) (Q ^ QU+ i 2k + 6 2 fe fc!r(fc + 6) 



^ 2fc 2'/! r(2fc + l + 6») > l " ' 2fc T(2fc + 1 

11=1 l = u x ' v 
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Let's denote p = u + 1, then we have m k — Yl P =x A k , P {9)6 p , where 
, ... 2k + 9 2 k klT(k + 9) . 



2k r(2k + l + 6) ' 
a ™ V- 1 2k + 6 2 k k\ T(k + l + 9) . 
i=p— i 



[PROOF OF LEMMA [2]: 

Proof. We use mathematical induction to show this. When p = 1, then 
2fc + 6> 2 fc fc! 2 fe - 1 (fc-l)! 



4^(0) 



2fc (2fc- l + 6»)---(fc + 6») (2fc- l + 0)---(fc + 0) 



and 



2 fe - 1 (fc - 1)! _ 2 k (k - ly.k! 
~ (2fc — 1) - - - A: ~ (2fc)! ~ M ' 



2 fc -Hfc-l)! 2*- 1 (fc-l)! _ 1 

Ml >~ (2fc-l + 0)-.-(fc + 0) ~ 2k---(k+l) ~ 2 M ' 



Moreover, Afe.i = ^^T^ < 1 < ^ = 2, therefore 

Now we assume that A k ^ p {9) satisfies the inequality 
^A ktP < A ktP {9) < A ktP < — 1 



(6) ^A k<p <A kiP (6)<A k>p <^,k>p>l. 



Since A k<p+1 (6) = J2i= P 2k W (2k-i+ey P -(k+i+e) ' then b y assumption @, we 
2p L> 2k 2H\ (2k - 1 + 1) • • • (k + I + 1) - fc ' p+1 1 J 



' ' 1 2 A /,-' .1/. 



- ^ 2k 2 l l\ (2k-l)---(k + l)' 



l, P 



Thus, 

1 ^ fc + ; 2 fc fc! r(fc + Q ^2*fe! r(fc + Q 

2 p2^ 2 fc 2 'Z! r(2fc + l) lp ~ A ww s 2^ 2 H\ r(2fe + l) 

1— P l—p ' 

Since ^ > h±k = i ; then 

1 ^fc + (2 fc )t! T(Hi) 1 k + l2 k kl T(k + l) 

2P^ 2k Wr(2fc + 1) '' p - 2^+1 2 fc 2'Z! T(2fc + 1) 



Hence, 

1 



2 p+i 



^4fc, P +i < ^4fe,p+i(^) < A kiP+1 . 
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We are going to show Ak, P < ^i^jVfc > p > 1. Since A kl i < ^P^-> tnen we assume 
that Ak.p < k > p, then 

_^ 2 fc fc! T(k + l) 1 ^ 2 k k!T(k + l) 

A k , P+ i - 2^ 2H\ fm + i) ,p - 2^2 ^ irfi+i)' 

Let's define = ^zfr^fc+i) ' fc — l>Z>p>2. In order to figure out the 
monotonicity of B k j in Z for fixed k. We should consider 

2 ; Z!r(fc + Z + l) fc + Z 



-B/tj 2 ; +!(Z + l)!r(fc + /) 2(2 + 1) " 

If Z < k - 2, then %^±i > 1, therefore we know B kk _ 2 OT B k k _i is the maximum 
term. Since 

2 fc fc! r(2fc-2) 1 

B k ,k-2 — 



2 k - 2 (k - 2)! r(2fc + 1) 2fc-l 
2 fe fc! r(2fc-l) 1 

i>fe,fe-l 



2 fe - 1 (/c- l)!T(2fc + l) 2fc — 1' 
then 

y- 1 2 fc fc!r(fc + Q fc-p _ fc-p 1 

^ 2 z nr(2A: + 1) ~ 2k - 1 ~~ 2(k-p) + 2p- 1 < 2' 

Thus, ^4fc, p +i < 2p^2^ = 2^t- F° r tnc ^ as ^ P art J we can also use mathematical 
induction to prove it! 

I a m A ,_, 2 k k\T(k + 0) 2 k k\T(k) 

lAk ' m AkM ~^2kv{2k + e) r(2fc + i)' 

2 k k\T(k) /T(2k)T(k + 6) 

~ T(2fc + i)lr(2fc + e)r(fc) ~ 

(2fc-l)-..fc 
S fc ' l| (2fc + 0-l).-.(fc + 0) 1 

= ^ 1 -2TT^T ) - (1 -^ ) - li 

2fe-l j ^2-1 

= MmI _ § *+* + fc <,<?< 2fe -i ft +«>)(*»+') 

+(-l) fc - 



(2fc - 1 + 6») • • • (fc + 0) 

->'• ■ -| #2-1 



+- 



(2ft - 1 + 6) ■ ■ ■ (k + 0) 

2fe-l 



.1Z2 (2fc — 1) - - - fc 1 

k<h<l 2 <2k-l 11 y > 

^A M |(i + ^L_)... ( i + I)-i| 
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- . 2/c 2/c — 1 k ~\~ 1 _ , „ . 

= 6A k J 1 = 6A kl 2 -l\=0A k i. 

' 1 2A; — 1 2fe — 2 k 1 1 1 

Therefore, \A k ,i(9) — Ak,x\ < 9A k ,i- Now we assume that 
(7) \A k , p (9)-A k J<0A ktP , 
then for \A ktP+1 (0) - A kjP+ i\, we have 

\Ak, P +i(6) — A k , p +i\ 

_^ 2k + 9 2 k k\T(k + l + 9) A ^ 2 fc fc!r(fc + Q A 

2k 2HT(2k + 1 + 6) lA j ^2H\T(2k + l) l ' pl 
i—p ?= p 

_^2k + 2 k kW(k + l + 6) 



2k 2 l l\T(2k+ 1 . 

i= P 

y2k + 8 2 k k\T(k + l + 9) y2 k k\T(k + l) 

+ l ^ 2k 2H\r(2k + l + 6) l ' p ~ ^ 2'Z!r(2fe+l) ' ,pl 

l—p l—p 

_ k ^2k + 9 2 k kir(k + i + e) 



fc-i 



2 fe fc!r(fc + o , r(jfe + / + 0)r(2fc) 



, 2^!r(2fc + 1) ,pl r(& + or(2fc + 1 

By the assumption ([7]), we have 

\Ak tP +i{0) — j4fe,p+i| 
! 1 2 fe fc!r(/c + I + 9) 
i= P 

k-l 



~ Up 2^ 2 l l\2kT(2k + 9) l ' p 



Moreover, 



2 fc fc!r(fc + 1) T(k + l + 8)T(2k) 
^ 2HW{2k+ 1) l ^T(2k + 9)T(k + l) 



r(fc + z + 6>)r(2fc) (2fc-i)---(fc + Q . 

T(2fc + 6>)r(fc + 1 ~ \2k + 6-l) ■■■[k + l + 9) 

-IC-j^T+s'-C-S+T+s)- 1 ! 



2fc-l 



E E 



+ --- + (-l) fe 1 

1 ; (2fc-l + 61) ■■■(k + l + 9) 1 
2k-i 1 

d \~ E 7^ + E (/ + 0wj 

l = k+l k+l<h<l 2 <2k-l y > y ' 

+ --- + (-l) fe 



(2fc-l + 6>)---(jfc + Z + . 
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2k — 1 

^ o gJ + --- + (2fc -i) 1 ..( fc + o l ^ l(1 + 2fc L T ) --- (1 + feT7 ) - 11 

_ fl | 2fc 2fc -! fc + ^ + ! 1| = g, 2fc _ g fe - ^ cg 

1 2fc — 1 2fc — 2 fc + Z 1 'fc + / 1 k + l 

and 

r(fc + z + 0) _ i 

2fcr(2fc + 6») _ 2fc(2fc + 0-l)---(fc-M + 0) 
1 T(k + l) 



< 



Therefore, 



2k{2k -!)■■■ {k + l) r(2fc + l)' 

2 k k\ T(k + l) 
2HI r(2fc+l)' 

2H\ T{2k + \) A ' 



, . , , . 2 fe fc! T(k + l) . 

\A k>p+1 - A kiP+1 \ < Op^ ol]] r( r )h _ uU A i,P 

l=p 

Thus, we have proved the lemma. 



[PROOF OF LEMMA [3]: 
Proof. By lemma [2j we have 



£ ,±<^il AkM < £ ,£.(">«*>" ' 



l=[\]+i fe=i J=[A]+i fc=; 
- 0^ » (Alog^ 

i=[A]+l fc=i 

, ^ i. Aloe . , ^ 0'~ A flW-^+i 2 

<4 £ (x)'e Alos ^=4 £ — =4 



2 ; ^2' 2W+ 1 2- 

i=[A]+l l=[X]+l 

0, as -> 0, due to [A] + 1 > A. 



[PROOF OF LEMMA [5]: 



□ 



□ 



Proof. By mathematical induction, we can prove this lemma. For p = 1, by Stir- 
ling's lemma 

V 2 e 

we have A/.^ ~ V^^f We can therefore assume that 

(8) A k , p ~ C p -£-(— ^—)'\ 

For Afc )P +i, we have 



_^ 2 fc fc!r(fc + Q 

felP+1 ~^ 2 ^!r(2fe + l) l ' p - 

l—p 
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By the assumption ©, we have Ve > 0, 3M > 0, such that Vfc > M, 
then rewrite Ak, P +i as X + Y", where 

2 fc fc!r(fc + Q 

^2*Z!r(2fc + l) ' ,p 

and 

^ 2 fc fc!r(fc + Q 
2^ 2 '/!r(2fc+i) 

/=M+1 v ' 

Define a fe (/) = 2' ;fr(2fc+i) ( ^+1 ) ' ' and Sl = £?=p yf °fe(0- Now we are g° in g 
show linife_j. +00 = 0, and lim^+oo = 1. Indeed 

max p <i< M (M -p+1) 2 k k\T(k + M) 



< X < 



2Pp\ T(2k + 1) 



maxp<i<M(M -p+1) j-k 
— — /it— 



M- 



2Pp\ v 4 fc ' 

and 

(1-e) J2 %k(l)<Y<(l + e) J2 

l=M+l l=M+l 



Then 



< X- < r x < ^ p <km(M + p + 1 , f 
" Si " ^a fc (ife) " 2Pp! V 4p j 



thus, lim/j^+oo = 0. Similarly, 



M 



c max p <i< M %(M-p + l) fc M-i 



then 



Ei 2Pp! v v 4p ; 

Therefore, 

d-Od-^^,^ < (1 + . )(1 _^^E,, 

ZjX Li 1 Lii 

letting k — >• +00, then e — > 0, we have 

y 

km — = 1, 

k— >oo 2-i\ 



hence A^ p ^i ~ Si. Moreover, 

2 fc fc!(fc-l)! 2(p+l) fc r(fc + Q p , _ p k 
k[) T(2fc + 1) [ p + 2 ' l\T(k) l 2(p+l) M 2(p+l) J ' 
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Let X^ be negative binomial, NB(k, a), where a = 2 (p+i) ' then 

_ 2^!(fc-l)! 2(p+l) fc k 

ak[l) r(2fc + i) ( « _ } ' 

and 

2*fcl(fc-l)l 2(p + l) fc ^ 1 fc _ 

r(2fc + i) ( ^T^ } 2^]f p ( x a-0- 

We claim that Y,tp -j p (X k = I) ~ lEt^K = 0, as fc ^ +oo, where 
Z = fca. Indeed, 



;7 



fc— 1 I — 

= E(\J l jrp(x k a = i) 



Ve > 0, consider 

fc-i 



E(Vf) Pp (*« = o 

£ {^7P{x k a = i)+ E {^7P{x k a = i) 



p<i<i (l-e) / (l-£)</<io(l+e) 



+ E (J l j) p p(x k = i), 



io(l+e)<'<fe-l 

then 



o< E {Myp{x k = i)<uj±y E p(x Q fc = o 
E (J\) p P{x k = i)<{J^Y E ^5 = 

-e)<Kfc-l V ^ J (l+e)<Kfc-l 



< 

i (l+e)<i<fc-l ' V r l (l+e)<l<k- 

= J l jyP(X*>l (l + e)) 



and 



) p P(loa-e)<X*<l o a + e)) 



lo{l + e) 



< E (^j) p p(x k = o 

lo(l-e)<«<J (l+e) 
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By LDP of NB(k, a), we have 



lirn (,fe"P(X*<io(l-e)) = 0, 



and 



Therefore, 



and 



lim (WV^>^o(l + e)) = 0. 

fe-i-+oo y p 



E (\jj) P P(X k a = -»• 0, as k -> +oo. 

p</<; (l-e) 



E (\[j) Pp (X k a = I) -> 0, as * 



+oo. 

V / 

Z (l+e)<Z<fc-l 

By central limit theorem of X^, letting fc — > +oo, then e — > 0, we have 
E (j!fY>P(XZ = l)^l, asfc^+oc. 

io(l-e)<Ki (l+e) 

Then 

1 4 fe fc!(fc-l)! /P + l, fc 
x ~ z | ° p r(2fc + l) l p + 2 j ' 

By Stirling's formula, we know 

( P + 2 y 1 ( £+j >|fc 



[PROOF OF LEMMAS]: 
Proof. Define 



s=0 

Then 



where follows binomial distribution B(k,f3), and /5 = ^=rr i 



A ' £ + 1 

Next, we show S 2 T^i^T + 7Tl) fe ' where s o = P k - To this end ' Ve > °> 



□ 



consider 

s=0 V 
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let's rewrite ((9]) as 



0<s<s (l-e) V (l-e)s <s<(l+e)s 



(l+e)s <s<fc-l 

Since 



and 



(l+e)s <s<fe-l 

then by LDP of we have 



sa. E <vW^=-'> 

0<s<s (l-e) 



and 



Moreover, 



fe-i>s>s (l+e) 



( '/ t - ( r + ^„ ) ' p((1 - e)8 °^^ (1 + e)a ° ) 



(l-e)«o<s<(l+<0«o 



E (a/^)^ = s ) 



< (J fc _ (1 _ e)so m(l - £ ) s o < Jf| < (1 + e)*o) 



By the central limit theorem of binomial distribution, 
lim P((l - 

k— >-\-oo 

Letting k —> +00, e — »■ 0, we have 



lim P((l - e)s <X%<(1 + e)a ) = 1. 

fc— >+oo ^ 



(l-£)s <s<(l+e)s 

Therefore, E 2 ~ (t^^tt^ + j^)*. Then we show C k ,i - E 2 . Since Ve > 
0, 3M > 0, such that Vfc - s > M, 

(k-s)2 W + 



where 



.sr a 

s=0 
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and 

E C)<V)— ■ 

k-M<s<k-l v 7 

Since 



and 



0<B< max A fe _ s ., ( " I (— —) k ~ M , 
k-M< s <k-i '\k-MJ 1+1 



B „ i . i / k \ , M 



Thus, linifc_j. +00 ^ = 0. Moreover, 



0< s <fe-M 
< A 



^ 1 + c > E H^'^tt^txt)*' 

^ V ° ' A (fc - s )5 I + 1 



0<s<fc-M 



and 



(l-e)(l- V ( fc W) S ^77^(7TT) & " S ) 



fc-M<s</c-I 
< A 



<(l + e)(l- E 0(^)^7^(7^)'-' 



k-M<s<k-l 



A y (fc-s)i v J + l' 



Since 



then 



~ E O^W^TTI^ 

k-M<s<k-l v 7 ^ */ 

< max 

- k-M<s<k-l L (fc-s)3 ^ + 1 J V s / A ' 

E( k\ ( \ — l\s /~< l_ \k— s 

k-M<s<k-i \s)\ x ) W+l 7 

q \ K ~ S ) 2 

— > as fc — » oo. 



Hence, letting — ► oo, e — ► 0, we have 



. A/ i 1 .A — I I , k 

= [ + (A-A)a + i) j2 fcI ( ^T + 7TT j 



□ 



[PROOF OF LEMMA [7]: 
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Proof. Let's assume that lime-j-o K n (9) — ( iL ^-) n , for 1) < A < u(u+l),u > 2, 
then we are going to show lim^o K n {9) = limg^o K n {9). K x {9) can be rewritten 
as 



where 



and 



Km = 



l + G x (6) ' 



Z^=i w Z^fc=z fci A k,i 

og 



Eaig , EZL t ^ jL (^.«(g)-^,i) 



We claim that lim^o = lim ^o G x {9) = 0, then 

IimiC„(0) = limKJd) 

is followed. Indeed, by lemma [21 we have 
< \F*(0)\ < 



Z)S=1 ^ Efctz (A '°fe| g) l^fc+w.iW _ A k+n,l\ 



< Z^=i w (W Z^fc=i Ei A 



fe+n,£ 



l^,i=i u Z^fe=i — m — Am 

Similarly, by lemma [2] we have 

0< |G A (6»)| < [A]0->O, as6»^0. 

Now we are going to show Let's assume that limg^o K n {9) = for u(u— 1) < 

A < u(tt + 1), it > 2. We can rewrite K n (9) as 

[A] /TO V^OO (A log ^) fc 4 r^ 00 ( A l°g /) 

E 17 2^ik=v fcl ^fc." 2^1 k=v fcl ^fc+W.D 

V^W fliv 00 ( A1 °S?) fc /I Y^°° (Alog^)^ ^ 

11=1 Z^/=l 17 Z^fc=i fe! Z^fe=-u fe! 

Then, 

Um K x {9) = V lim ^ k=v M- * M lim ^ k=v M f fc+ "'" . 

6(^0 ■^e->0v^[ A ] fli Y^°° ( A1 °g?) 4 61—5-0 v^oo (A log 7 ) fc ^ 

" =1 Z^i = l u l^k=l fe! Z^fe=„ fe! A k,v 

By lemma |4] and lemma [5j we know 

Y^oo (Alog^) fc , 
r lm Z^fc^-u fe! S±k+n,v _ , V . n 

0_>O Y"°° ( AlQ g?)" /I W + l ' 

l^k=v fe! fe,U 



then we need to show 

i^o W A 1 Y^°° x - ' s 

Z^l 17 2^fe=i fe! ^fe 



/it, v^oo (A log |) . 
i- 17 Z^fc=u fe! c / \ 
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Once we have this, then lim e ^ k n (6) = £ [ X i 1 8 (u _ 1} (v)(^) n = (^)". To 
end, we divide 9 X to the numerator and denumerator of 

(Alog£)* 



7 2^,k=v k\ ' 



thus, we need consider 



/ro-A v^oo (A log . 
l^k=v fcl ^fc," 

V [A] fli-AV™ ( A1 °s?)* 4, ,' 

Z^/ = l 17 2^fc = Z fcl 

Since 1 < w < [A] < A, then 



9 — w — = (z) — w — 



fc! v 6» 7 ^ fc! 

s— A:— ij 

y (Alog|)V 
2^ fc] °M» 



where 



Then 



' fc\ . A — v . 



s=0 



p«-A v^°° ( Xl °sj) k a sr°° ( A1 °s ?) fc n 

7 Z^fc=« fcl ^-fc.f _ l^k=v k\ ^ k ' v 



V*W fll-JV 00 (A log I)'' ^ ^[A] ^p 00 (Alog^)^ 

To figure out the limit, we must find the leading term among 

fc=; 

By lemma |4] and lemma [6j we have 



Y^oo (A logjQ f-j r , 

r 2^ik=v fci °M !■ U M 
llm oi M fc = llm ~FT- 

0^0 V^oo (A log ^ fc^+oo Ckl 

L~ik=l k\ ^ k > 1 



f-l j A Ii^^"l7 1 A-u i v 

y _ (\-v)(v+Tj ) l — + — , fc 

To find the leading term, we need to figure our the maximum term among ^j- 
j^, 1 < I < [A]. Consider /(a:) = ^ + ^ - 2 - (f + ^), then 

1 {X) = (z + l) 2 ~ A' 

we know 

f>0, ie<VA-1 
\< 0, a; > a/A - 1. 
Therefore, + -p^ attains its maximum at [vA] — 1 or [vA]. 
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Case 1: For (u — l)u < A < u 2 , [y/X] = u — 1, since 

/( „- 2) = 2 -(^ + -l T ) 

and 

/(«-l) = 2-(^i + I), 
A u 



then f(u — 2) — /(w — 1) = i — m , < 0, so f(u — 1) is the maximum term. 



j,(u-l) 
X-l i j 
r A , . ' +1 . < 1,V1 < Z < [Al. 

then 

lim 7 ^-=5„_i(0,Vl<i<[A]. 

Case 2: For u 2 < A < u(u + 1), [\/A] = u, then the maximum term of + ^ 
should be /(it — 1) or f(u). Since 

■^- 1 )-^4-^TT)> ' 

then f(u— 1) is the maximum term, and 

, , - 1-1 /• A. 



a-; i j_ 

a z+i 



Thus, 

lim 7 ^-=5„_i(0,Vl<i< [A]. 

Case 3: For A = u(u + 1), [\/A] = w, then the maximum term of + ^ i 
f(u— 1) or /(u). Since 



is 



/(„-!)-/(„) = I- ! 



then 



and 



But 



a-; _| i_ 

A <l,Vl</< [A],Z^u,u-l, 



A— u+l i u_ 
A 



A— U+l I M— 1 

A tt _ 2 

A— ^ I u 

A " r u+l 



,. Cfe, u (1 + (A-«)(u+l)) 2 1 n 

hm — — — = hm - irrr— = 0, 

^+°° c m-i (i+ (A ^; 1 1)M )— fc^ 

therefore, Ck, u -i is the maximum term. We have 

lim 7 ^-=<5 u _i(0,l</<[A]. 

-to y^oo (A log ^) fc . 

l im ^k= v fc! AM = j (u) 1 < -y < [A] . 



Thus, 



Z^/=l 17 Z^fc=Z fc| 71 



□ 
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